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FOREWORD 


This  report  is  based  on  theoretical  studies  presented  in  Naval  Surface  Warfare 
Center  (NAVSWC)  Acoustics  Group  Letter  Reports  LR-109  entitled  "Effective 
Dynamic  Material  Parameters  of  Inhomogeneous  Materials.  I.  Preliminary  Results 
Based  on  the  T-Matrix  for  a  Pair  of  Spherical  Scatterers,”  of  21  March  1989  and  LR 
117  entitled  "Effective  Dynamic  Material  Parameters  of  Inhomogeneous  Materials. 
II.  Multiple  Scattering  and  Mode  Conversion,”  of  10  September  1990. 

Effects  of  localized  multiple  scattering  on  the  longitudinal  wave  vector  of 
coherent  waves  in  inhomogeneous  matei  ials  are  studied.  A  new  wave  vector  is 
derived  by  replacing  the  single  particle  far  field  forward  scattering  amplitude  by  the 
average  pair  scattering  amplitude.  The  Perkus-Yevick  hard  sphere  radial 
distribution  function  is  used  to  average  over  the  separations  of  scatterers.  The 
results  are  applied  to  mass  loaded  polyurethane  spheres  in  syntactic  foam,  a  type  of 
material  potentially  significant  in  Naval  applications  and  to  lead  spheres  in  epoxy 
resin.  In  the  case  of  lead  in  epoxy,  it  is  found  that  the  pair  scattering  mechanism 
significantly  alters  the  real  and  imaginary  parts  of  the  effective  wave  vector  and  the 
mode  conversion  scattering  cross  section.  Only  negligible  effects  are  observed  in 
mass  loaded  polyurethane  in  syntactic  foam. 

This  work  was  sponsored  by  the  NAVSWC  IR  Program  and  the  Mines  Block 
Program/Mine  Supporting  Technology /Covert  Materials  (Code  U077). 

The  author  wishes  to  thank  Mrs.  Linda  Tsui  (E52)  for  transferring  our  'Pair 
Scattering  Algorithm’  to  the  Naval  Research  Laboratory  Supercomputer  and 
running  it  to  obtain  the  data  documented  in  this  report. 
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INTRODUCTION 


The  present  paper  is  based  on  NAVSWC  Acoustics  Group  Letter  Reports  LR- 
109  entitled  "Effective  Dynamic  Material  Parameters  of  Inhomogeneous  Materials. 

I.  Preliminary  Results  based  on  the  T-Matrix  for  a  Pair  of  Spherical  Scatterers” 
dated  3/21/89  and  LR-117  entitled  "Effective  Dynamic  Material  Parameters  of 
Inhomogeneous  Materials.  II.  Multiple  Scattering  and  Mode  Conversion”  dated 
9/10/90.  The  essential  elements  of  this  work  were  also  reported  at  professional 
society  meetings.1 

As  stated  in  LR-109,  the  purpose  of  the  work  was  to  explore  the  possibility  of 
including  in  the  classical  theory  for  the  effective  wave  vector  certain  localized 
multiple  scattering  events  which  have  thus  far  been  ignored.  Replacing  the  single 
particle  longitudinal  far  field  forward  scattering  amplitude,  fp(l,0),  by  the 
corresponding  (averaged)  pair  scattering  amplitude,  fp(2,0),  achieves  this  objective. 
Effects  of  mutual  dynamic  loading  of  scatterers,  as  well  as  multiple  scattering 
between  nearest  neighbors  will  be  accounted  for  and,  if  important,  will  show  up  in  the 
effective  parameters,  particularly  as  the  concentration  of  scatterers  is  increased.  P 
should  however  be  realized  that  the  classical  theory  strictly  applies  only  to  low 
concentrations  of  inclusions  and/or  weak  scattering.  In  the  present  work,  we  are 
mainly  interested  in  high  concentrations  and  strong  scattering  such  that  the  product 
of  the  number  density  of  inclusions,  n,  times  the  extinction  cross  section,  oe,  is  much 
larger  than  unity  near  resonance.  This  means  that  power  loss  per  unit  length  due  to 
scattering  and  absorption  is  much  greater  than  1/e,  the  coherent  wave  is  strongly 
attenuated  and  the  classical  theory  not  strictly  applicable.  The  present  results  must 
therefore  be  compared  with  experiment  to  test  their  validity.  This  cannot  be 
emphasized  strongly  enough.  Using  the  far  field  scattering  amplitude  in  a  dense 
medium  and  neglecting  possible  near  field  effects,  as  is  done  here,  although 
conceptually  attractive,  can  at  best  only  be  a  heuristic  approach  to  the  problem. 

In  LR-109  we  presented  an  analysis  of  the  algorithm  for  the  T-Matrix  of  a  pair 
of  scatterers,  the  latter  being  the  essential  ingredient  in  this  type  of  calculation. 
Preliminary  estimates  of  the  effective  attenuation  and  wave  speeds  of  longitudinal 
waves  in  inhomogeneous  syntactic  foam  with  silicone  inclusions  and  in  rubber  with 
air  bubbles  were  also  presented.  The  calculations  were  preliminary  in  the  sense  that 
only  the  spatial  orientation  average  of  a  pair  of  scatterers  was  accounted  for  in  fp(2,0) 
and  constant  nearest  neighbor  distances  were  assumed  for  given  volume 
concentrations,  the  latter  being  representative  of  a  regular  lattice  rather  than  a 
random  medium.  There  was  also  a  problem  with  the  convergence  of  the  algorithm 
whenever  the  maximum  mode  number  exceeded  certain  low  values.  Only  very 
restricted  ka-ranges  could  be  explored. 

The  convergence  problem  (due  to  the  large  integer  storage  limitation  of  the 
NAVSWC  CDC  860  computer  (see  also  LR-117))  was  solved  subsequently  and  the 
algorithm  applied  to  two  different  inhomogeneous  systems;  mass  loaded 
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polyurethane  in  syntactic  foam  and  lead  in  epoxy  resin  (NAVSWC  Acoustics  Group 
Letter  Report  LR-117).  The  new  ingredient  in  the  latter  work  was  the  use  of  the 
Perkus-Yevick  pair  distribution  function  for  the  statistical  distribution  of  nearest 
neighbor  separations  in  a  random  distribution  of  spherical  inclusions  in  an 
inhomogeneous  solid. 

The  body  of  the  present  paper  is  essentially  a  reprint  of  NAVSWC  Acoustics 
Group  Letter  Report  LR-117.  The  Appendix  was  taken  from  NAVSWC  Acoustics 
Group  Letter  Report  LR-109.  The  other  data  presented  in  LR-109,  although  correct, 
is  incomplete  because  of  the  convergence  problem  cited  above  and  has  therefore  not 
been  included  in  the  present  summary. 
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THE  MODIFIED  CLASSICAL  EXPRESSION 
FOR  THE  EFFECTIVE  WAVE  VECTOR 


The  classical  expression  for  the  effective  wave  vector  of  longitudinal  waves  in 
inhomogeneous  media  is:2 


k*  =  k  +  2nnf(l,0)/k  (1) 

p  p  p  p 

where  kp  =  co/cp,  cp  is  the  longitudinal  wave  speed  of  the  matrix  material  and  fp(  1,0  ) 
the  far  field  longitudinal  forward  scattering  amplitude  of  a  single  scatterer. 
Replacing  n  by  the  volume  concentration  of  spherical  scatterers  of  radius  a, 

<$>  =  4nna3/3,  we  have: 


k*  =  k  [1  +  (3<b/k2a3)f  d,0)/2]  (2’ 

P  P  P  P 

* 

Equations  1  and  2  show  that  Akp  s  (kp  -  kp)  is  proportional  to  n  times  the  far 
field  scattering  contribution  of  a  single  scatterer  scattering  independently.  We  wish 
to  replace  this  term  by  a  quantity  proportional  to  fp(2,0,<l>),  the  volume  concentration 
dependent  average  of  the  far  field  forward  scattering  amplitude  of  a  pair  of  scatterers 
where  the  average  has  to  be  performed  with  respect  to  all  orientations  and  over  a 
range  of  separations  of  scatterers.  In  so  doing,  the  modified  version  of  the  equation 
will  account  the  mutual  dynamic  loading  of  scattered  in  close  proximity  and  for 
all  orders  of  mutual  scattering  events  between  pairs. 

The  orientation  average  of  fp(2,0),  the  far  field  forward  scattering  amplitude  of 
an  isolated  pair  of  scatterers,  may  be  performed  analytically.  The  result  is  given 
below.  For  the  radial  average  we  use  the  volume  concentration  dependent  Percus- 
Yevick  hard  sphere  radial  distribution  function,  g(4>,r).3  Accordingly,  the  average 
number  of  scatterers  with  centers  within  dr  at  the  radial  distance  r  from  any  other 
scatterer  is 


N(<t»,r)  =  4nng<4>,r)r'dr 

Since  we  only  wish  to  consider  unobstructed  pairs  of  scatterers  in  fp(2,0),  those 
not  having  intervening  neighbors,  we  have  to  choose  a  cutoff  radius  for  g(4>,r).  In  the 
present  work  we  have  explored  the  first  nearest  neighbor  maximum  of  g(4>,r),  Rq,  as  a 
cutoff.  Defining  the  probability  of  finding  the  center  of  another  scatterer  within  dr  at 
the  radial  distance  r  from  any  other  scatterer  as 


The  number  of  centers  in  the  shell  between  rand  r  •+  dr 

P(4>.r»dr  =  - 

The  number  of  centers  in  the  shell  between  2a  and  R 
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we  have 

gi^.rir'dr  g(4),x)x'dx 

Fl4>,r)dr  =  — -  =  — -  =  pi^.xidx  1 3 ' 

10  ■)  ,  /  *0  •> 
g  <  <Jj  ,rlr“d  r  g<4>,x)x"d\ 

•I  2a  I  I 

where  x  =  r/2a  and  xo  =  Ro/2a.  The  average  pair  scattering  amplitude  per  scatterer 
at  the  volume  concentration  4>  is  therefore 


f  (2,0,4))  —  [  p)4)A)  f  (2,0,x)dx  (4) 

P  j  l  P 

where  fp(2,0,x)  is  the  average  of  fp(2,0)  over  all  directions  of  the  incident  wave  with 
respect  to  the  axis  of  the  pair.  The  average  pair  scattering  amplitude  per  unit  volume 
at  the  volume  concentration  is  therefore 

F  (2,0,4))  -  (n/2)f  (2,0,4>)  u5 1 

p  p 

where  the  factor  1/2  accounts  for  the  fact  that  on  making  the  change  from  the  average 
pair  scattering  amplitude  per  scatterer  to  that  per  unit  volume,  all  pairs  are  counted 
twice. 

Replacing  the  quantity  nfp(  1 ,0)  in  Equation  ( 1 )  by  Fp(2,0,4>)  we  obtain  the 
modified  classical  wave  vector 


k  =  k  +  nnf  {2,0,4>)/k  (6 

p  p  p  p 

or 


f  (  f  (2,0,x)  g(4),x)x2dx 
—  •  ri  J  1  P 

k  =  k  +  (34>/4k  a J)  -  (7) 

P  P  P  x 

r  o  •> , 

J  g(4>,x)xax 

which  shows  explicitly  that  the  linear  dependence  of  the  effective  wave  vector  on  in 
Equation  (2)  has  now  become  nonlinear  in  terms  of  g(<J>,x).  From  Equations  (2)  and 
(7),  we  obtain  expressions  for  the  reduced  wave  speed  and  the  attenuation  coefficient 
at  given  volume  concentrations.  In  general,  if  k  is  the  wave  vector  of  the  matrix  and 
k  that  of  the  composite,  we  have 

k  =  (uv  c  I  +  i  a 

(  c/c)  =  k/Reul  (  k  )  (8) 

and 


Attenuation/unit  length  =  8  69  Im(  k  )  -  -lOlogiuk’  21ni1  k  )  (9) 
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where,  in  Equation  (8),  we  have  assumed  that  the  matrix  medium  is  elastic.  Plots  of 
Equations  (8)  and  (9)  are  given  below  for  single  particle  and  pair  scattering  at 
different  volume  concentrations  of  inclusions. 
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LONGITUDINAL  FAR  FIELD  FORWARD  SCATTERING 
AMPLITUDES  AND  SCATTERING  CROSS  SECTIONS 


The  longitudinal  forward  scattering  amplitude  of  isolated  spherical  scatterers 
in  terms  of  the  T-Matrix  is.4,5 


f  (1,0)  =  (1/ik  )  V  (2n  +  l)TU  (10 

p  P  ~  <' 

n  =0 

where  Tn11  is  the  element  Tnn  11  of T",  a  diagonal  matrix  (see  the  Appendix  for  more 
details).  In  addition  to  the  attenuation  (extinction)  coefficient  and  the  wave  speed 
which  follow  from  Equations  (1),  (8),  (9),  and  (10),  it  is  also  of  interest  to  study  the 
absorption  coefficient  and  the  scattering  cross  sections.  There  are  three  types  of 
scattering  cross  sections  which  are  of  interest  here,  opp,  the  longitudinal  to 
longitudinal  wave  scattering  cross  section,  o^s,  the  longitudinal  to  shear  wave  (mode 
conversion)  cross  section  and  oa,  the  absorption  cross  section.  For  isolated  spherical 
scatterers  the  first  two  of  these  are4,5 


o  (1  )  =  (4ii/k2) 
PP  p 


<2n+  1 ) 

n  =  0 


111) 


O 


(1)  =  (4u/k2)  V  (2n  +  1 ) 

pt,  p  — 

n  =  0 


(12) 


The  absorption  cross  section  is  non-zero  only  if  the  inclusion  material  is  viscoelastic. 
We  assume  the  matrix  material  to  be  elastic.  The  absorption  cross  section  may  be 
obtained  by  means  of  the  Forward  Scattering  Theorem,  according  to  which  the 
forward  scattering  amplitude,  fp(kp,kp)  is  related  to  the  cross  sections  as4  6,7 

A  A 

(4n/k  )Im|f  (k  ,k  |  =  o  +o  +o=o  (13) 

p  ppp  pppsae 

where  oe  is  the  extinction  (attenuation)  cross  section.  Equation  (13)  is  perfectly 
general  result,  applying  to  any  matrix/scatterer  combination  and  scatterer  shape, 
including  any  cluster  of  scatterers,  such  as  a  pair  and  regardless  of  whether  the 
scatter  is  elastic  or  viscoelastic.  In  fact  it  applies  to  any  kind  of  scattering 
phenomenon  and  is  used  extensively  in  such  diverse  fields  as  astronomy,  nuclear  and 
elementary  particle  scattering,  radar  propagation,  medical  diagnostics,  etc.  Using 
Equations  (10),  (11),  and  (12)  in  Equation  (13)  we  obtain  oa(l)  for  single  scatterers. 

The  absorption  coefficient  follows  from  the  definition 


-  no  1 1 i 

Absorptiorvunit  length  =  —  10 log  .J e  “  )  =  4  343no  (1)  (14) 

10  H 


/ 


NAVSWCTR  91-182 


With  the  aid  of  Equations  (1),  (9),  and  (13)  this  may  be  written 


Absorption/unit  length  =  AttenuatioiVunit  length -4. 343n[o  (l)  +  o  ( 1 ) j  (15) 

pp  pt> 

for  single  scatterers.  For  pairs  of  scatterers  we  simply  replace  n  by  n/2  and  the  single 
particle  scattering  cross  sections  oPn(l),  OpS(l),  oa(l)  by  the  corresponding  pair 
scattering  cross  sections  opp(2),  ops(2)  and  oa(2). 

Since  Equation  (13)  holds  for  any  orientation  and  any  separation  of  a  pair  of 
scatterers  it  also  holds  for  the  averaged  quantities 

(4n/k  )Im(f  (2 ,0 ,<b )  1  =  o  +  o  +  o  (16) 

P  p  PP  Pb  a 

where  the  bars  denote  orientation  plus  separation  averages.  Hence  the  average 
absorption  cross  section  for  a  pair  follows  from  Equation  (16)  and  the  average 
absorption  coefficient  from  Equation  (14)  with  n  replaced  by  n/2  and  oa(l)  replaced  by 
oa.  In  our  plots  of  cross  sections  we  have  normalized  those  for  single  particles  by  na2 
and  those  for  two  particle  averages  by  2na2.  The  latter  does  not  represent  the  true 
average  of  the  geometrical  cross  section  of  a  pair  but  the  definition  is  convenient  since 
the  two  cross  sections  should  (and  do)  approach  each  other  at  very  large  particle 
separations. 

As  mentioned  above,  the  Perkus-Yevick  average  over  radial  separations  must 
be  done  numerically,  while  that  over  angular  orientations  of  the  axis  of  the  pair  with 
respect  to  the  incident  plain  wave  direction  may  be  done  analytically.  The  result  of 
the  orientation  average,  assuming  a  uniform  angular  distribution,  in  terms  of  the 
T-Matrix  elements  of  a  pair  of  scatterers  is 


f  (2,0,x)  =  il/ik  )  lim 
P  P, 

max-** 


j 

/  max 

(  X 


m  =0 


! 

max  -  m 
V 

1  =  0 


ni 


rr.ee  (111) 

ml,!/ 


(17) 


O 

PP 


(4  ii/k2) 
P 


I  i 


/  max 

max  -  m 

.  ,  O  . 

(  I 

\  e 

— -  m 

|  ml,!  |  / 

m  =  0 

i.r  =  u 

(18) 


o  =  (4n/k2)  lim  ( 

/  max 

max  -  m 

V  c 

1  |  rpOfc*  (41  ) 
III  -  , , 

1  2 

1  + 

Tee<5n  |2j 

pf>  P  ,  ' 

m 

[  !  ml,!’ 

1 

m  1  ,  1 

1 

max— 

m  =  0 

i,r=o 

'  J 

(19) 


where  em  =  2-8m,o-  The  analysis  of  the  T-Matrix  of  a  pair  of  scatterers  in  terms  of  the 
T-Matrices  of  the  individual  scatterers  is  given  in  Reference_8.  Substituting 
Equation  (17)  into  Equations  (4)  and  (7)  yields  ^(2,0,4>)  and  kp*,  respectively. 
Averaging  oPp  and  ops  in  Equations  (18)  and  (19)  with  respect  to  p($,x)  from  x  =  1  to 
x  =  xq  yields  tne  averaged  cross  sections  opp  and  ops  in  Equation  (16). 


The  super  and  subscripts  on  the  matrix  elements  are  interpreted  as  follows:  In 
the  Appendix  we  show  that  the  T-Matrix  has  the  form 
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T 


XU 

0 

0 

x14 

XI5 

0 

0_ 

*22 

x23 

0 

0 

x26 

0 

x32 

x33 

0 

0 

x36 

X41 

0 

0 

x44 

x45 

0 

X51 

0 

0 

x54 

x55 

0 

0 

x62 

x63 

0 

0 

x66. 

(20) 


where  the  xjj  are  'm-diagonal’  matrices  (Equation  (4)  of  the  Appendix).  The  2x2  block 
in  the  upper  left  hand  corner  corresponds  to  longitudinal  scattering,  the  4x4  block  in 
the  lower  right  hand  corner  to  shear  wave  scattering,  all  others  to  mode  conversion 
scattering.  In  particular 


v  _  rpeeiHi 
II  1 


_  fptie  (41 1 
X41  "  1 


_  met  ( 51  ( 

51  - 

The  added  superscripts  on  the  T-Matrix  refer  to  the  parity  of  the  corresponding  wave 
functions,  e.g.,  T°e<4*’  transforms  an  even  parity  longitudinal  wave  function  into  an 
odd  parity  shear  wave  function.  .The  first  sum  in  Equations  (17),  (18),  and  (19)  is  the 
sum  over  the  azimuthal  index  or  the  'm-blocks’.  Subsequent  sums  are  over  the 
elements  in  the  'm-blocks'  or  the  mode  numbers.  lmax  is  the  largest  mode  number 
included  in  the  calculations.  In  the  present  work  lmax  =  9,  where  1  =  0  corresponds  to 
the  monopole,  1  =  1  to  the  dipole,  1  =  2  to  the  quadrupole,  etc. 

In  the  case  of  elastic  materials,  the  T-Matrix  must  satisfy  the  equation 
(T-Matrix  Theorem)4 


0  =  Realm  +  tV  (21) 

where  the  superscripts  t  and  *  denote  the  transpose  and  the  complex  conjugate 
respectively.  Since  T  is  symmetric  in  the  case  of  elastic  materials  this  equation  may 
also  be  written 


0  =  Real(T)  +  TT* 

Doing  the  multiplication  indicated  in  Equation  (21)  shows  that  the  diagonal  elements 
of -Real  [Tee(U)],  the  elements  in  Equation  (17),  are,  in  fact,  equal  to  the  combined 
sums  over  T  of  the  elements  in  Equations  ( i  8)  and  (19),  for  fixed  m  and  1.  This 
illustrates  the  fact  that  the  Forward  Scattering  Theorem  for  elastic  scatterers  is  a 
corrolary  of  the  T-Matrix  Theorem.  In  the  case  of  single  spherical  scatterers,  the 
situation  is  somewhat  simpler  algebraically,  as  may  be  seen  from  Equations  (12)  and 
(13)  in  the  Appendix. 
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In  the  case  of  viscoelastic  inclusions  Equations  (17),  (18),  and  (19)  are  still  true, 
butT1  *  T  and  the  left  hand  side  of  Equation  (21)  is  not  equal  to  zero.  As  a  conse¬ 
quence,  the  absorption  cross  sections  in  Equations  (13)  and  (16)  are  not  equal  to  zero. 

Since  Equations  (10),  (11),  (12),  (17),  (18),  and  (19)  are  true  regardless  of 
whether  the  inclusions  are  elastic  or  viscoelastics,  the  same  algorithm  was  used  in 
both  cases  to  calculate  the  various  cross  sections,  the  attenuation  coefficient,  the 
wave  speed,  and  the  absorption  coefficient  in  the  case  of  viscoelastic  scatterers.  In  the 
case  of  elastic  inclusions  in  elastic  matrix  materials,  we  have  checked  the 
convergence  of  our  pair  scattering  algorithm  by  plotting  Equation  (9),  using  the  right 
and  left  hand  sides  of  Equation  (16),  as  defined  by  Equations  (17),  (18),  and  (19),  with 
6a  set  to  zero.  The  two  curves  should  be  (and  are)  identical.  This  shows  that  in  the 
case  of  lead  in  epoxy,  choosing  lmax  =  9  is  sufficient  for  the  calculation  of  the 
T-matrix  (strictly  only  for  its  real  part)  in  the  range  from  kpa  =  .1  to  1.  Figure  1 
shows  the  error.  Comparing  with  the  corresponding  attenuation  curves  shows  that 
the  error  is  of  order  10  4  at  worst  (at  K1*A  =  1.  (=  kpa)  for  $  =  5.2%).  In  the  case  of 
viscoelastic  inclusions,  the  difference,  as  calculated  for  a  viscoelastic  polymer  in 
syntactic  foam,  is  due  to  absorption. 
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FIGURE  1  THE  DIFFERENCE  BETWEEN  THE  ATTENUATION  COEFFICIENTS 
(IN  dB/m)  DERIVED  FROM  THE  RIGHT  AND  LEFT  HAND  SIDES  OF 
THE  T-MATRIX  EQUATION  FOR  LEAD  IN  EPOXY  (EPON  828Z) 
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RESULTS  AND  DISCUSSION 


PRELIMINARIES 

To  begin  with,  in  Figure  2  we  have  plotted  the  three  Perkus-Yevick  distribu¬ 
tions  used  in  the  present  work.3  Note  the  relatively  large  oscillations  at  high  concen¬ 
trations  of  inclusions  and  the  large  numbe*  of  scatterers  in  contact  at  57.6  percent 
compared  to  the  number  at  5.2  percent.  These  differences  determine  the  nonlinear 
dependence  of  kp  on  $. 

All  calculations  are  for  scatterer  radii,  a(  *  A)  of  2  mm.  All  plots,  except  those  of 
the  radial  distribution  function,  are  versus  Kl*A,  where  K1  ( =  kp)  is  the  wave  vector 
of  the  respective  matrix  material.  The  effective  wave  speeds  of  inhomogeneous 
media,  <C1  >,  have  been  reduced  by  those  of  the  corresponding  matrix  media, 

Cl(-cp). 

The  material  parameters  of  the  components  of  the  two  inhomogeneous 
materials  discussed  in  this  work  are  listed  in  Table  1.  All  parameters  are  assumed  to 
be  constant.  The  polyurethane  based  elastomer  should  be  thought  of  as  being  mass 
loaded  with  metal  dust  in  order  to  achieve  the  mass  density  given.  It  should  also  be 
thought  of  as  containing  about  10  percent  by  volume  of  microscopically  dispersed  air 
for  enhanced  absorption.  The  order  of  magnitude  of  the  shear  loss  tangent  of  the 
polyurethane  inclusions  is  typical  of  the  shear  loss  tangent  peak  in  the  relaxation 
region  of  this  type  of  modified  elastomer. 


TABLE  1.  MATERIAL  PARAMETERS 


p(gm/cm3) 

A* 

Go* 

8* 

Syntactic  Foam 

0.64 

0.938 

0.344 

0  0 

Polyurethane 

2.64 

0.320 

0.015 

0.45 

Epoxy  (Epon  828Z) 

1.20 

2.185 

0.769 

0.0 

I^ead 

11.3 

17.10 

3.714 

0.0 

*  Moduli  are  reduced  by  2. 25E10  dyne/cm2 
G  =  G0(l-i8)  =  Shear  Modulus 
A  =  First  Lame  Parameter 


The  frequency  corresponding  to  K1*A  =  .1  in  the  syntactic  foam  matrix 
(sound  speed  =  2391  m/sec)  is  19  kHz  (7.6  kHz  for  A  =  .5  cm),  that  in  the  epoxy 
matrix  (sound  speed  =  2640  m/sec)  is  21  kHz  (8.40  kHz  for  A  =  .5  cm). 
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FIGURE  2  THE  "HARD  SPHERE,"  PERKUS-YEVICK  RADIAL  DISTRIBUTION 
FUNCTION,  g<4>,x)J  X0(<|>)  =  Ro(<t>)/2A  IS  THE  "CUTOFF”  RADIUS 
USED  IN  THE  PRESENT  CALCULATIONS 
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CROSS  SECTIONS 

The  salient  features  of  the  frequency  dependence  of  the  effective  wave  speed, 
attenuation  and  absorption  may  be  understood  best  by  first  considering  the  various 
cross  sections. 

Figures  3  and  4  show  that  as  the  volume  concentration  increases,  the  pair 
scattering  cross  sections  broaden  in  the  low  frequency,  dipole  (mode  conversion) 
region  and  the  dipole  maxima  decrease.  The  effect  is  most  conspicuous  in  the  case  of 
Pb  in  epoxy,  Figure  3,  but  appears  to  saturate  above  31.4  percent.  Similar  remarks 
apply  to  urethane  in  syntactic  foam,  where  the  low  frequency  maximum  is  due  to  the 
dipole,  the  high  frequency  maximum  to  the  monopole  and  intermediate  maximum  is 
due  to  the  quadrupole,  Figure  4.  The  absorption  cross  section  is  a  significant  part  of 
the  total  in  the  latter  case,  Figure  4  and  the  difference  between  opp  +  ops  and  opp 
shows  that  mode  conversion  is  significant  in  both  the  dipole  and  £ne  quadrupole 
regions  of  scattering.  In  the  case  of  5.2  percent  Pb  in  epoxy,  the  dipole  maximum  is 
clearly  shifted  to  slightly  higher  frequencies.  At  higher  concentrations  the  shift  is 
insignificant. 

Physically,  the  decrease  of  the  scattering  cross  section,  ops,  may  be  thought  of  as 
due  to  the  changing  relative  phase  of  motion  of  the  two  parts  of  the  pair  of  scatterers 
as  the  frequency  of  oscillation  changes.  As  fare  as  the  dipole  resonance  is  concerned, 
the  pair  is  simple  a  double  oscillator,  constrained  by  the  matrix  material  and  coupled 
by  multiple  scattering.1  In  a  coupled,  linear  oscillator  the  single  particle  resonance 
splits  into  two;  at  the  lower  resonance  the  two  components  tend  to  move  in  phase,  at 
the  higher  resonance  they  tend  to  move  out  of  phase  with  respect  to  each  other. 
Depending  on  the  orientation  of  the  axis  of  the  two  particles  with  respect  to  the 
incoming  wave,  these  resonances  are  excited  more  or  less  strongly.1-^  Destructive 
interference  and  partial  cancellation  of  the  shear  waves  scattered  by  each  of  the  two 
scatterers  takes  place  in  the  resonance  region,  resulting  in  a  net  decrease  of  the  total 
scattered  shear  wave  amplitude.  Below  the  resonances  the  two  parts  of  the  oscillator 
tend  to  move  in  phase  with  respect  to  the  matrix  material.  Above  the  resonances 
they  move  out  of  phase  with  respect  to  the  matrix  material. 

The  average  monopole  and  quadrupole  oscillations  of  the  two  parts  of  the 
oscillator  in  the  case  of  urethane  in  syntactic  foam  seem  to  occur  independently;  the 
scattering  cross  sections  are  largely  unefTected  by  considering  the  two  scatterers  as  a 
double  oscillator,  Figure  4,  although  in  the  case  of  the  monopole  our  limited 
resolution  does  not  allow  a  firm  conclusion.  Only  the  average  relative  phase  of  the 
center  of  mass  motions  of  the  two  parts  of  the  oscillator,  their  dipole  motion,  appears 
to  become  coupled  in  the  above  sense. 

We  know  from  other,  unpublished  studies  of  ours  that  the  detailed  structure  of 
the  dipole  resonance  region  depends  on  the  cutoff,  Xo.  The  major  effects,  however,  the 
broadening  and  the  decrease  of  the  maximum  of  the  mode  conversion  cross  section, 
are  invariant  with  respect  to  Xo.  As  Xo  decreases,  the  interference  minimum  induced 
by  the  high  frequency  out  of  phase  dipole  resonance  shows  up  clearly  at  high  volume 
concentrations  of  Pb  in  epoxy,  not  only  in  the  cross  sections  and  the  attenuation,  but 
also  in  the  wave  speed. 
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FIGURE  4  SCATTERING  CROSS  SECTIONS  OF  POLYURETHANE  SPHERES  IN 
SYNTACTIC  FOAM  ( — )  :  SINGLE  PARTICLES.  (-  O  A -) : 

PAIRS  OF  PARTICLES  AVERAGED  OVER  ORIENTATIONS  AND 
CENTER  TO-CENTERSPACINGS 
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WAVE  SPEED 

The  most  significant  effects  are  observed  in  the  case  of  Pb  in  epoxy,  Figure  5. 
The  dispersion  region  appears  to  broaden  slightly  in  the  pair  scattering  mode  with 
increasing  concentrations  of  inclusions  and  the  single  particle  maximum  between 
Kl*A  =  .4  and  .5  is  greatly  reduced.  These  effects  appear  to  be  due  to  the  phase  of 
relative  motion  as  discussed  above.  In  the  case  of  urethane  in  syntactic  foam  on  the 
other  hand,  multiple  scattering,  or  coupling  effects  appear  to  be  negligible  at  all 
concentrations  of  inclusions,  Figure  6,  again  as  expected  from  the  discussion  of  the 
cross  sections. 

Some  measurements  of  the  wave  speed  for  lead  in  epoxy  resin  are  available.10 
They  show  that  the  dispersion  minimum  in  the  wave  speed  shifts  to  higher 
frequencies  as  the  volume  concentration  of  inclusions  increases.  The  present  results 
do  not  show  this  effect. 


ATTENUATION 

At  low  frequencies,  attenuation  due  to  pair  scattering  is  enhanced  by  about  a 
factor  of  two.  Near  the  peak  of  the  dipole  resonance  it  is  decreased  somewhat, 
Figures  7  and  8,  as  expected  from  the  discussion  of  the  cross  sections.  The  effects  are 
most  conspicuous  for  Pb  in  epoxy.  Since  attenuation  and  absorption  are  proportional 
to  the  appropriate  cross  sections,  Figures  7,  8,  and  9  would  look  exactly  like  the 
corresponding  cross  section  plots  if  the  ordinates  were  linear.  Figure  3  therefore 
shows  clearly  that,  e.g.,  at  57.6  percent  Pb  in  epoxy  peak  attenuation  in  the  pair 
scattering  mode  is  reduced  by  about  25  percent.  The  difference  between  the  two 
curves  for  pair  scattering  in  Figure  8  is  due  to  absorption,  as  explained  above. 


ABSORPTION 

Comparing  Figures  8  and  9  shows  clearly  that  absorption  in  the  viscoelastic 
polyurethane  inclusions  is  a  substantial  part  of  the  total  attenuation  mechanism  in 
this  materia!  at  all  concentrations  of  inclusions;  absorption  and  scattering  make 
approximately  equal  contributions.  Except  at  very  low  frequencies,  either  way  of 
calculating  absorption  leads  to  essentially  the  same  result,  Figure  9.  The  irregular 
behavior  of  the  points  based  on  pair  scattering  in  the  low  frequency  region  is  not 
understood  at  present.  It  may  be  due  to  computational  inaccuracies  resulting  from 
the  finite  shear  loss  tangent. 
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FIGURE  5.  THE  EFFECTIVE  WAVE  SPEED  IN  EPOXY  (EPON  828Z)  CONTAINING  LEAD 
SPHERES,  CALCULATED  ON  THE  BASIS  OF  SINGLE  PARTICLE 
SCATTERING  ( — )  AND  PAIR  SCATTERING  AVERAGED  OVER 
ORIENTATIONS  AND  CENTER-TO-CENTER  SPACINGS  (-  □  -) 
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FIGURE  6  THE  EFFECTIVE  WAVE  SPEED  IN  SYNTACTIC  FOAM  CONTAINING  POLY 
URETHANE  SPHERES,  CALCULATED  ON  THE  BASIS  OF  SINGLE 
PARTICLE  SCATTERING  ( — )  AN  1)  PAIR  SCATTERING  AVERAGED 
OVER  ORIENTATIONS  AND  CENTER  TO-CENTER  SPAClNGSl  □  -) 
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FIGURE  7. 
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THE  EFFECTIVE  ATTENUATION  IN  EPOXY  (EPON  828Z)  CONTAINING  LEAD 
SPHERES,  CALCULATED  ON  THE  BASIS  OF  SINGLE  PARTICLE  SCATTERING 
( — )  AND  PAIR  SCATTERING  AVERAGED  OVER  ORIENTATIONS  AND 
CENTER-TO-CENTER  SPACINGS  f- □  )  (-  0  -) :  ATTENUATION  FOR 
AVERAGED  PAIR  SCATTERING  PREDICTED  BY  opp  +  opi> 
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FIGURE  8.  THE  EFFECTIVE  ATTENUATION  IN  SYNTACTIC  FOAM  CONTAINING  POLY¬ 
URETHANE  SPHERES,  CALCULATED  ON  THE  BASIS  OF  SINGLE  PARTICLE 
SCATTERING  ( — )  AND  PAIR  SCATTERING  AVERAGED  OVER  ORIENTATIONS 
AND  CENTER  TO  CENTER  SPACINGSI- □ -).  (-0-):  ATTENUATION  WITHOUT 
ABSORPTION  FOR  AVERAGED  PAIR  SCATTERING 
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CONCLUSIONS 


Whereas  some  of  the  amplitude  changes  predicted  by  using  pair  scattering  in 
place  of  single  particle  scattering  in  the  case  of  Pb  in  epoxy  are  large  enough  to  be 
measurable,  in  the  case  of  the  polyurethane-syntactic  foam  system  the  effects  are 
probably  in  the  noise,  except  perhaps  for  the  level  of  low  frequency  attenuation. 
Moreover,  no  significant  shifts  in  the  frequencies  of  wave  speed  dispersion  and  the 
attenuation  edge  are  predicted,  except  in  the  case  of  Pb  in  epoxy  where  the 
attenuation  edge  moves  to  somewhat  lower  frequencies. 

The  general  absence  of  significant  frequency  shifts  in  the  observable  quantities 
points  toward  the  fact  that  the  splitting  of  the  dipole  resonance  frequency  due  to 
coupling  between  the  two  oscillators  of  the  pair  of  scatterers  is  weak.  In  the  case  of 
Pb  in  epoxy,  however,  this  splitting  appears  to  contribute  toward  moving  the 
attenuation  edge  to  somewhat  lower  frequencies.  Destructive  interference  effects  due 
to  the  high  frequency  resonance  of  the  double  oscillator  lead  to  generally  decreased 
amplitudes  of  observable  quantities  in  the  resonance  region.  Both  frequency 
splitting  and  destructive  interference  are  responsible  for  the  broadening  of  the  dipole 
resonance  region,  particularly  in  the  case  of  Pb  in  epoxy. 
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APPENDIX 

AN  ALGORITHM  FOR  THE  T-MATRIX  OF  A  PAIR  OF 
SPHERICAL  SCATTERERS  OF  ELASTIC  WAVES 


The  T-matrix  to  be  calculated  is:A  1 

Tp  =  Rl(d)T(l)(l  -o(2d)T(2)ot(2d)T(]  *^R(d)  +  o(2d)T(2)Rl(d)^  + 

R(d )T(2 )  ( 1  - o l(2d  )T(  1  )o(2d )T(2) )  " '  (  Rl(d )  +  ol(2d )T( ]  )R(d ) )  ( A- 1 ) 

where  2d  is  the  distance  between  the  scatterers,  T  is  the  single  particle  T-matrix  and 
R  and  o  are  translation  operators. 

We  have  analyzed  this  matrix  in  some  detail  and  have  used  the  insights  gained 
to  maximize  the  computational  efficiency  of  our  algorithm. 

Of  special  interest  is  the  case  T(1)  =  T(2)  for  spherical  scatterers  with  relative 
translations  along  the  z-axis.  In  that  case,  we  show  that  Equation  (A-l)  takes  the 
form 


T  =  A  +  PAP  (A-2) 

p 

where  A  is  the  first  term  in  Equation  (A-l)  and  PAP  the  second  (or  vice  versa).  P  is 
given  below.  Since  P  is  a  very  simple  matrix,  essentially  only  one-half  of  Tp  needs  to 
be  calculated  and  only  one  inversion  needs  to  be  performed. 

The  analysis  was  carried  out  in  the  space  in  which  the  polarization  and  parity 
indices  are  not  mixed.  Our  algorithm  also  runs  primarily  in  this  space.  We  call  this 
Space  I.  For  computational  efficiency,  inversions  must,  however,  be  carried  out  in  the 
space  in  which  these  indices  are  mixed. A  2  We  call  the  latter  space  Space  II.  The  two 
spaces  are  related  by  a  similarity  transformation. 

Using  Pao’s  wave  function  convention/  3  the  R  and  o  matrices  for  translations 
along  the  z-axis  have  the  following  form: 
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where  each  xy  represents  an  "m-diagonal”  system  of  square  matrices,  the  one  in  the 
upper  left  hand  corner  being  the  largest  (of  dimension  lmax  + 1,  where  Imax  is  the 
maximum  mode  index)  corresponding  to  the  azimuthal  index  m  =  0,  and  the  one  in 
the  lower  right  hand  corner  being  a  scalar  (corresponding  to  the  azimuthal  index 
lmax)' 


Every  block  of  x,  x(m),  has  the  form 
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where  sjj  and  ay  are  complex  or  real  numbers  (for  o  and  R,  respectively)  and  they  have 
the  property  that  in  the  transpose  of  x(m),  xl(m),  ay  has  changed  sign  but  sy  has  not: 
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where  ay  =  -aji  and  sy  =  Sjj. 

Hence,  for  any  xy  in  Equation  (A-3) 


xl  =  I)-1  x  D 
>j  ij 

where  D  has  the  form 


(A-7) 
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D  = 


(A-8) 


and  D2  =  I  or  D  1  =  D.  This  matrix  may  be  thought  of  as  a  string  of  diagonal 
"m-blocks”  along  the  main  diagonal  with  alternating  plus  and  minus  ones,  each  block 
corresponding  to  some  x(m)  and  starting  with  + 1  in  the  upper  left  hand  corner. 


The  structures  of  RT,  oT,  R'T,  olT,  oT(2)olT(l),  otT(l)oT(2),  etc.,  in  fact  the 
structures  of  all  of  the  matrices  in  the  various  products  in  Equation  (A-l),  including 
that  of  Tp  itself,  are  all  the  same  and  are  of  the  form: 
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where  every  x„-.matrix  has  the  form  of  Equation  (A-4)  and  where  in  oT  and  RT, 
X23  =  X14  =  0.  Similarly  for  o'T  and  RlT.  Focussing  on  oT(2)olT(  1 ),  one  finds  that  a 
typical  xjj  has  the  form: 


N 

x  =  Y  o(ik)T(2 jk)otUk)T(l,nk) 
k  =  1 

where  1<N<3  and,  e.g.,  o(ik)  is  an  element  of  the  set 


I  «ee 
(°n  ’ 


ee  _  00  _  eo  ' 
°”  1  °22  ’  °22  ’  °23  J 


11 


(A-10) 


(A-l  1) 


and 


_ee  _ 
°H=X11 


.00  =. 


,eo  =. 


,0,1  X22  ’  °22  X  33  X55’n22  X  44  X  66’°23  X  36  X4.c  X54 


=  —  x 


63 


in  Equation  (A-3)  and  all  of  these  have  the  form  of  Equation  (A-4).  Also  T^  is  an 
element  of  the  set 


(T  T  T  T  T  ) 
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and  it  should  be  noted  that  T  has  the  form: 
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(A-13) 


where  T]  1  =  xn  =  X22,T22  =  x33  =  x44,T33  =  X55  =  X6b,Ti3  =  X15  =  X26,T3i  = 
X51  —  x62  in  Equation  (A-13)  and  all  of  the  Tn  are  diagonal  (for  spheres).  In  the  case 
of  elastic  scatterers,  Tjj  =  Tjj1,  i.e.,  T13  =  T31  tor  elastic  spheres. 

For  olT(l)oT(2)  we  have  matrices  of  the  form: 


N 

x  =  ±  ^  ot(ik)T(ljk)onk)T(2,nk) 

k  =  1 


(A-14) 


The  sign  in  Equation  (A-14)  is  positive  if  there  are  an  even  number  of  o23eu  in  every 
term  and  negative  otherwise.  It  can  be  shown  that  if  oT(2)olT(  1 )  has  the  form  shown 
in  Equation  (A-9),  then  olT(l)oT(2)  has  the  form: 
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where  the  sings  are  those  in  Equation  (A-14).  The  exact  relationship  between  a 
typical  xn  in  Equation  (A-9)  and  a  typical  xjj  in  Equation  (A-15)  may  now  be  found  as 
follows.  Using  Equation  (A-7): 


o(ik)T(2Jk)olGk)TU  ,nk)  =  DDoliJ  Ol)T(2  I)o(lk>  DT(lfnk)I)D  = 

Dol(ik)T(2jk)o(lk)T(l,nk)D  (A-16) 

since  the  T-matrices  are  diagonal  (for  spheres). 

Hence, 

x  =  ±  D  x  D  and  x  =  ±  D  x  D  (A - 1 7 ) 

‘j  ij  ij  ij 
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if  T(l)  =  T(2).  Therefore,  if  T(l)  =  T(2),  by  assembling  the  D  matrices  into  a 
"supermatrix,”  P,  which  also  handles  the  sign  changes  on  going  from  Equation  (A-9) 
to  Equation  (A- 15),  we  find  that: 


P 


o 


\ 


Pi 


(A- 1 8) 


and 

(1  -o*ToT)  =  P(1  -oTo*T)P  (A-19) 

Similarly, 

(l-o^'oTT1  =  P(1 -oTotT)_,P  (A-20) 

i.e.,  (l-aTolT) 1  transforms  exactly  like  (l-oTo^).  Similar  considerations  show  (by 
inserting  P2  between  the  matrices)  that  P  transforms  any  part  of  the  first  term  of 
Equation  (A-l),  i.e.,  any  pat  of  A,  including  A  itself,  into  the  corresponding  part  of  the 
second  term  and  vice  versa.  Hence,  in  analogy  with  Equation  (A-20),  the  second  part 
ofTpis: 


B  =  PAP 


(A-21 ) 


and 


Tp  =  A  +  B  (A-22) 

if  T(l)  =  T(2).  This  result  depends  only  on  the  spherical  geometry  of  the  scatterers 
[Equation  (A-16)],  not  on  their  elastic  properties.  The  spheres  may  be  elastic  or 
viscoelastic. 


SYMMETRIES  AND  COMPUTATIONAL  SAVINGS 

Upon  studying  the  forms  of  oT(2)olT(l)  and  1  )aT(2)  more  closely,  one  finds 
with  respect  to  the  matrix  elements  in  Equation  (A-9),  for  m  <  lmax  + 1 


*22  =  X11  ,  *32  - X41  ,  x62  ~  X51 


x14 - x23  •  x44  -  x33  ,  x54 - x63 


(A-23) 


x26  -  x15  >  x36  - x45  .  x66  -  x55 

and  similarly  for  xjj  in  Equation  (A-17)  and,  in  fact,  for  all  of  the  matrices  and  matrix 
products  in  Equation  (A-l),  including  Tp  itself.  When  m  =  lmax+ 1  one  finds 
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x22  ~  x32  -  x62  -  x  1 4  -  x44  -  x54  -  x26  -  x36  -  x66  = 


x4 1  =  x63  ~  x23  =  x45  =  0.  (A-24) 

except  for  T(lj  and  T(2),  where  Equation  (A-23)  holds  for  all  m. 

If  one  transforms  from  Space  I  to  Space  II,  where  in  every  m-diagonal  sub-block 
we  count  through  all  polarization  and  parity  indices  at  fixed  mode  number,  1 
(forming  a  6*6  block  for  each  1 )  and  repeat  that  count  for  all  1  <  m  (in  every 
m-diagonal  block),  we  find  that  the  above  symmetries  hold  between  adjacent  columns 
of  the  transformed  matrix.  This  is  the  form  of  the  matrix  (l-oTo'T)  that  should  be 
inverted  because  it  is  m-diagonal  and  may  be  inverted  one  block  at  a  time,  yielding 
much  smaller  matrices  than  in  Space  I.  But  because  of  the  above  symmetries 
between  the  elements  it  is  not  necessary  to  invert  the  whole  matrix.  Only  every  other 
column  needs  to  be  inverted;  i.e.,  only  50  percent  of  every  m-diagonal  block  in 
Space  II  needs  to  be  inverted  and,  as  shown  above,  only  one  inversion  is  necessary 
when  T(l)  =T(2),  reducing  the  computational  effort  by  another  50  percent  for  a  total 
reduction  of  computations  associated  with  matrix  inversion  of  75  percent  when 
T(1)  =  T(2)  and  50  percent  when  T(1)*T(2). 

Since  the  symmetries  indicated  in  Equations  (A-23)  and  (A-24)  hold  throughout 
all  matrix  manipulations,  once  RT,  oT,  RlT,  and  o^  have  been  formed  it  is  only 
necessary  to  work  with  the  nine  blocks  (1,1),  (4,1),  (5,1),  (2,3),  (3,3),  (6,3),  (1,5),  (4,5), 
and  (5,5)  in  Space  I.  The  other  nine  non-zero  blocks  may  be  ignored.  Since  there  are 
36  blocks  altogether,  we  thus  only  work  with  one  quarter  of  the  matrices,  for  a 
considerable  savings  in  computational  effort  (T(l)  =  T(  2)  or  not).  These  nine  blocks 
correspond  to  the  odd  numbered  columns  of  (l-oTolT)  in  Space  II,  which  are  the  ones 
being  inverted. 

It  should  also  be  noted  that  in  Space  I  the  m-diagonal  blocks,  x(m),  in  Equa¬ 
tion  (A-4)  are  densely  populated  with  non-zero  elements,  whereas  in  Space  II  all  of 
the  zero  blocks  in  Equation  (A-9)  are  carried  and  stored.  It  is  therefore  preferable  to 
work  in  Space  I  when  matrix  multiplications  are  performed  since  many  zero  manipu¬ 
lations  are  avoided  and  only  the  indicated  subset  of  nine  matrix  blocks  is  considered. 
This  then  is  the  space  in  which  most  of  our  calculations  are  being  performed.  On  the 
other  hand,  for  matrix  inversions,  Space  II  is  indispensible. 


THE  STRUCTURE  OF  THE  ALGORITHM 

1.  Do  the  matrix  multiplications  to  arrive  at  (1-oTo'T)  in  Space  I. 

2.  For  the  inversions  of  (1-oTo'T)  do  the  similarity  transformation  which  takes 
this  matrix  into  Space  13.  This  amounts  to  nothing  more  than  a  repacking  of  the 
matrix  into  m-diagonal  blocks  and  requires  no  greater  computational  effort  than 
loading  it  in  the  first  place. 

3.  Invert  very  other  column  of  this  matrix.  If  T(1)*T(2)  do  the  same  for  (1-oToT). 

4.  Do  the  inverse  similarity  transformation  to  take  (1-oTo'T) 1  back  into  Space  I. 
Again,  this  amounts  to  no  more  than  an  unpacking  of  the  matrix  into  the  nine 
working  blocks  of  Space  I.  IfT(l)*T(2)  do  the  same  for  (l-olToT)  ’. 
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5.  Perform  the  remaining  matrix  manipulations  in  Space  I.  For  T(1)  =  T(2), 
calculate  either  the  first  or  the  second  part  of  Equation  (A-l),  then  use 
Equation  (A-21)  to  obtain  Tp. 

By  properly  sequencing  the  computations,  storage  space  may  be  minimized. 
Also,  in  the  case  of  elastic  spherical  scatterers,  it  is  only  necessary  to  calculate 

R'T(l)  (1  -oT(2)olT(l))_1oT(2)Rl  (A-25) 

or 

RT(2)  (1  -olT(l)oT(2))_1olT(l)R  (A-26) 

since  one  is  the  transpose  of  the  other. 

The  results  may  be  checked  by  observing  that: 

1.  Tp  must  be  symmetric  for  elastic  materials  (true  for  any  T-matrix  when 
absorption  is  zero). 

2.  The  T-matrix  theorem  must  hold,  i.e., 

-ReaKT  )  =  T  T*  (A-27) 

p  p  P 

3.  In  the  special  case  ofT(  l)  =  T(2)  we  note  from  Equation  (A-2)  that 

T  =  PT  P  (A-28) 

p  p 

where  Tp  has  the  structure  of  Equation  (A-9),  is  m-block  diagonal  in  every  xij  and 
every  Xy  has  the  form  of  Equation  (A-5).  Hence,  since  this  transformation  changes 
the  sign  of  every  other  off-diagonal  element  in  every  m-block,  Equation  (A-28) 
implies  that  every  other  off-diagonal  element  in  very  m-block  of  Tp  must  be  zero. 

It  should  also  be  noted  that  the  transformation  from  A  to  B  in  Equation  (A-21) 
or  from  (l-oTolT)  to  (l-oToT),  etc.,  reduces  to  the  nine,  simple  manipulations: 

ill  =  Dxi)l),  X41  =  -DX41I)  ,  X5i  =  DX51D 
*23  =  -DX23I)  ,  X33  =  DX33D  ,  Xg3  =  -DX63D 
*15  =  DX15D  ,  X45  =  0x45!)  ,  X55  =  DX55D  . 
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